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Abstract 

A significant proportion of patients with multiple sclerosis who receive interferon beta (IFNP) therapy develop neutralizing 
antibodies (NAbs) that reduce drug efficacy. To investigate if HLA class I and II alleles are associated with development of 
NAbs against IFNP we analyzed whether NAb status and development of NAb titers high enough to be biologically relevant 
(>150 tenfold reduction units/ml) correlated with the HLA allele group carriage in a cohort of 903 Swedish patients with 
multiple sclerosis treated with either intramuscular IFNp-1 a, subcutaneous IFNp-1 a or subcutaneous IFN p-1 b. Carriage of 
HLA-DRB1*15 was associated with increased risk of developing NAbs and high NAb titers. After stratification based on type 
of IFNP preparation, HLA-DRB1*15 carriage was observed to increase the risk of developing NAbs as well as high NAb titers 
against both subcutaneous and intramuscular IFNp-1 a. Furthermore, in patients receiving subcutaneous IFNp-la carriage of 
HLA-DQA1*05 decreased the risk for high NAb titers. In IFNp-1 b treated patients, HLA-DRB1*04 increased the risk of 
developing high NAb titers, and in a subgroup analysis of DRB1*04 alleles the risk for NAbs was increased in DRB1*04:01 
carriers. In conclusion, there is a preparation-specific genetically determined risk to develop NAbs against IFNP high enough 
to be clinically relevant in treatment decisions for patients with multiple sclerosis if confirmed in future studies. However, 
choice of IFNP preparation still remains the single most significant determinant for the risk of developing NAbs. 
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Introduction 

Recombinant interferon beta (IFNP) is a first-line therapy for 
relapsing-remitting multiple sclerosis (MS) and its therapeutic 
effect includes reduced frequency of clinical relapses and new 
lesions [1,2] and reduced accumulation of disability over time [3]. 
Two types of recombinant IFNP products are used for treatment 
of MS, IFNP-la and IFNP-lb. The IFNP-la preparations are 
injected either at 30 u.g intramuscularly (i.m.) once a week or at 
22 (0.g and 44 p.g subcutaneously (s.c) three times a week, and the 
IFNP- lb preparation is injected at 250 U.g s.c. every other day. 
IFNP-la is identical in amino acid sequence to human IFNP, 



while IFNP-lb lacks the N-terminal methionine (Metj) and the 
asparagine-linked glycosylation at position 80 (Asnso) and has a 
cysteine-to-serine substitution at position 17 (CySi 7 Ser). 

A proportion of MS patients receiving IFNP treatment develop 
neutralizing anti-drug antibodies (NAbs), which at high titers block 
the biological response [4,5] and the therapeutic efficacy [6-8] of 
IFNp. The frequency of NAb positive patients differs depending 
on the IFNP preparations used with IFNP-lb inducing NAbs in 
28%-47% of treated patients, s.c. IFNP-la in 5%-28% and i.m. 
IFNP-la in 2%-6% [9]. Among the NAb positive patients, the 
proportion with high NAb titers has been reported to be 19%- 
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73% if treated with IFNfi-lb, 52%-79% if treated with s.c. IFNfJ- 
la and 19%-79% if treated with i.m. IFNfi-la [10-15]. 

Differences in immunogenicity between preparations are likely 
due to factors such as dosage, injection frequency, application site, 
glycosylation and protein aggregate formation, and particle 
contents [9,16,17]. Since not all IFNfi-treated patients develop 
NAbs, other patient-specific immunological factors seem to 
influence this risk. The human leukocyte antigen (HLA) region 
on chromosome 6p2 1 carries the genes with the strongest effect on 
MS susceptibility and has also been linked to the immunogenicity 
of IFNfi therapy. An association was shown between the 
DRB1*15:01-DQB1*06:02 haplotype and increased proliferative 
response of CD4+ T-cells to specific IFNfi epitopes in vitro [18]. 
Furthermore, a study on peripheral blood mononuclear cells from 
IFNfi-treated MS patients found an association between the HLA- 
DRB1*07:01-DQA1*02:01 haplotype and IFNfS immunogenicity 
[19]. Recent studies based on larger patient materials have 
associated an increased risk of developing anti-drug antibodies 
against IFNfi with the HLA-DRB1*04:01 and HLA-DRB1*04:08 
alleles [20,21]. Also, genetic associations of an intergenic single 
nucleotide polymorphism (SNP) on chromosome 8q24.3 
(rs4961252) and a SNP within the HLA region (rs9272105), 
independent of DRB 1*04:01 and DRB1*04:08, with the develop- 
ment of high anti-drug antibody titers against IFNfi have been 
reported [22]. 

The aim of this study was to investigate if HLA class I and II 
allele groups are associated with risk of developing NAbs in a 
Swedish cohort of IFNfi-treated MS patients and if there is an 
HLA influence on biologically relevant high NAb titers. 

Materials and Methods 

Ethics statement 

The study obtained approval from the ethical board of 
Karolinska Institutet, Stockholm, Sweden, and subsequently 
genotyped blood samples were obtained only from patients who 
gave their oral and written informed consent in accordance with 
the ethical permissions. 

Patients included in the study 

In June 2012, the Swedish NAb biobank of routine serum 
samples contained more than 7000 samples from more than 4000 
unique patients. Cross-referencing with the genetic database 
resulted in 1 726 patients with both NAb status and HLA genotype 
data, either by classical typing (HLA-A, HLA-C and HLA-DRB1 by 
Olerup SSP™ HLA Low resolution Kit [23] and HLA-B by 
Luminex based reverse SSO Labtype® SSO from One Lambda 
Inc., Canoga Park, CA, USA) or imputed either by using 
HLA*IMP:01 [24] with genotypes from the IMSGC WTCCC2 
MS genome wide association study (GWAS) [25], or with 
HLA*IMP:02 [26] using genotypes from the Immunochip [27]. 
The threshold of the quality score Q_ for using imputed genotypes 
was 0.7. In cases where several sources of a HLA genotype were 
found, classical genotyping data was primarily used, and for 
imputed data Immunochip was predominantly used whereas 
GWAS derived data was only occasionally used. Comparison of 
genotypes between these sources yielded 1068 pairs of which 952 
(89.1%) were concordant for both alleles, whereas 68 (6.4%) were 
discordant for one or both alleles. For 48 of the 1068 (4.5%) pairs, 
one allele was concordant but the other allele remained unknown 
from imputation due to the quality cut off. Genotypes for HLA- 
DQA1 and HLA-DQB1 were derived only from imputation of SNP 
data (Table SI). Classically sub typed genotypes and imputed 
genotypes from GWAS and Immunochip SNP data were used to 



identify HLA-DRB1*04:01 and HLA-DRB1*04:04 positive pa- 
tients. NAb positive patients were defined as those ever having a 
positive sample, and patients ever having a NAb titer above 150 
tenfold reduction units per milliliter (TRU/ml) were classified as 
high positive. NAb negative patients were required to never being 
tested positive and having at least one negative sample taken after 
twelve months of treatment. After exclusion of patients of non- 
Scandinavian origin (n = 403), defined as self-reported or by 
principal component analysis (PCA) performed within other 
projects [25,27]and patients not fulfilling the criteria for NAb 
positivity and negativity (n = 420), 364 NAb positive and 539 NAb 
negative patients were included in the study. Information on the 
type of IFNfi preparation used was available for 358 of the 364 
NAb positive patients and for all NAb negative patients (Table 1). 

Analysis for neutralizing anti-drug antibodies 

The antibody neutralizing activity in serum was assessed at the 
NAb laboratory at Karolinska Institutet in Stockholm. Patient 
samples sent to the laboratory in 2003-2006 were analyzed with 
the myxovirus resistance protein A (MxA) assay (MPA) [15] while 
later samples were analyzed with the MxA gene expression assay 
(MGA) as previously described [14,28]. A strong correlation of 
NAb titers was found when 1 1 8 serum samples were titrated with 
the MPA and the MGA assay in our laboratory, and the methods 
can thus be considered to yield similar results (Spearman r = 0.93, 
95% confidence interval 0.90 to 0.95, p value<0.0001). In both 
assays, titers were adjusted using the Kawade formula [29] in 
Softmax Pro (Softmax Pro Software, Sunnyvale, CA) and 
expressed as TRU/ml. Patients' NAb status was classified 
according to the following titer categories: negative (<10 TRU/ 
ml), positive (>10 TRU/ml) and high positive (>150 TRU/ml). 
High NAb titers > 1 50 TRU/ml markedly affect the in vivo 
bioactivity of IFNfi therapy by reducing the MxA gene expression 
in treated patients [5]. Titers > 150 TRU/ml were therefore 
considered as biologically relevant titers. 

Treatment groups 

Treatment information was obtained from the Swedish Multiple 
Sclerosis registry (SMS-reg; http://www.msreg.net) with permis- 
sion from the board. Patients were categorized based on the IFNfi 
preparation (i.m. IFNfi-la, s.c. IFNfi-la or IFNfi-lb) used at the 
time of defined NAb status. NAb positive patients and those with 
biologically relevant titers were categorized based on the 



Table 1. Characteristics of the cohort. 



Total NAb positive * NAb negative 



Number 903 364 539 

Age (mean years, IQR) 48.4(39-57) 51.2 (43-58) 46.5 (38-55) 

Female (%) 69.9 65.9 72.5 

Titers > 150 TRU/ml 209 209 0 

IFNp-1a i.m. 346 49 297 

IFNp-las.c. 355 172 183 

IFN|!-1b 196 137 59 



a IFNp preparation could be assigned to 358 of the 364 NAb positive patients, 

while the type of preparation used was unspecified for 6 NAb positive patients. 

Abbreviations: IQR = interquartile range, i.m. = intramuscular, 

s.c. = subcutaneous, TRU/ml -tenfold reduction units per milliliter, 

IFN[} = interferon beta, NAb = neutralizing antibodies 

doi:1 0.1 371 /joumal.pone.0090479.t001 
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preparation used at the sampling date of the first positive sample 
and the first sample with a titer > 1 50 TRU/ml, respectively. 

Allele frequency and carriage of alleles 

The frequencies of all analyzed HLA alleles are presented in 
Table S2 and S3. Within each treatment group, all allele groups 
with a count of at least ten carriers in the positive and negative 
group were analyzed for association with development of NAb as 
well as with biologically relevant titers. 

Analysis of HLA-DRB1*04:01 and HLA-DRB1*04:04 

Since not all DRB1 *04 positive patients were subtyped to 4-digit 
resolution, the association analysis of specific DRB1*()4 alleles to 
NAb development had to be weighted to compensate for this lack 
of information. A univariate association analysis was done by 
comparing carriers of a certain DRB1*04 allele against the 
DRB1*04 negatives together with DRB1*04 positives who were 
subtyped but who did not carry the allele in question. Hence, 
DRB1*04 positive individuals who were not subtyped were 
excluded from analysis. Since high resolution data were available 
only for a subset (88.4%) of the DRB1*04 positive patients and the 
proportion of subtyped patients was uneven between the NAb 
positive and NAb negative group (Table S6), this difference 
influences the result of a normal association analysis. To overcome 
this problem the reference population (DRB1*04 negative) was 
weighted when calculating the association of DRB1*04 alleles. 
This weighing constant is based on the percentage of subtyped 
individuals within the NAb positive and NAb negative groups and 
keeps the proportions of DRB1*04 positive intact through the 
analysis. For example, for the NAb positive individuals, 105 out of 
120 DRB1*04 positive individuals were subtyped (87.5%), and this 
group was thus compared with 87.5% of the number of DRB1*04 
negatives (in this case 207.4 of 237). 

DRB 1*04:08 could not be imputed due to low frequency in the 
reference panel used by HLA-IMP:01/IMP:02. The Swedish 
genetic material from the Epidemiologic Investigation of Rheu- 
matoid Arthritis (EIRA) study [30], which also has been geno typed 
with the Immunochip and imputed with HLA*IMP:02 and have 
classical DRB1*04 subtyping data, was used as a cohort for 
checking HLA imputation accuracy. We found that 33 out of 36 
known DRB1*04:08 positive individuals were imputed as 
DRB1*04:01 and three individuals were imputed as 
DRB1*04:04. Since both DRB1*04:01 and DRB1*04:08 have 
been identified as risk alleles for NAb development [21] it was 
decided to use and analyze the imputed data as it was, even 
though the estimate of DRB 1*04:01 will be too high and 
DRB1*04:08 too low. Among the 68 classically genotyped patients 
there was only one patient carrying DRB1*04:08 (1.5%). Of the 
185 patients with imputed high resolution data for DRB 1*04 we 
estimate that two or three patients (based on the percentage found 
in classical typing) are assigned DRB1*04:01 instead. Since we 
know that the frequency of DRB1*04:08 is incorrectly estimated 
due to the problems with imputation, and the frequency obtained 
from classical genotyping was low, we decided not to analyze an 
association to DRB1*04:08 in detail. The single classically typed 
DRB1*04:08 positive patient was NAb negative and received s.c. 
IFNP-la. 

Statistics 

All calculations were done in R version 2.15.1 [31]. Fisher's 
exact test (fisher.test in R) was used for all comparisons; two-tailed 
P-values were reported and considered significant if below alpha 
0.05. The Bonferroni method was used to correct P-values for 
multiple comparisons (Pc). 



Treatment specific absolute risk of NAb development 

The absolute risk (AR) for development of NAbs and 
biologically relevant titers was used to assess the influence of 
HLA allele group carriage within each of the three treatment 
groups. The AR was estimated using Bayes' theorem [32] : 



P(N + \G)-- 



P{N+)P{G\N+) 



P(N+)P(G\N+)+P(N-)P(G\N- 



where G was the allele group carriage and N + and N~ represented 
NAb development or negativity (or presence of biologically 
relevant titers and negativity). P(N) was the frequency of NAb 
development or biologically relevant titers within each treatment 
group regardless of genotype. 

HLA binding prediction to IFN(3 

The NetMHCII and NetMHCpan prediction servers [33,34] 
were used for binding prediction of recombinant IFNfS-la (166 
amino acids) and recombinant IFN(3-lb (165 amino acids) to HLA 
class II alleles DRB1*04:01, DRB1*04:04 and DRB1* 15:01, and 
the HLA class I alleles B*07:01 and C*07:01. The IFNP sequences 
were obtained from DrugBank (www.drugbank.ca; Accession 
numbers: DB00060 and DB00068). 

Results 

Association of HLA alleles with NAb development 

An association to DRB1*15 was observed for both development 
of NAbs (OR 1.43; P c = 0.036) and biologicaUy relevant NAb titers 
(OR 1.58; P„ = 0.01). Carriage of HLA-B*07 was associated with 
development of biologically relevant titers (OR 1.61; P c = 0.0497) 
(Table 2). Because of the high linkage disequilibrium (LD) within 
the HLA, a stratification on carriage of DRB1*15 was performed 
to assess whether B*07 was independendy associated to the 
development of biologically relevant titers. Since B*07 did not 
show an independent association within either stratum we assume 
that this association is due to LD with DRB1*15. All analyzed 
HLA allele groups can be found in Table S2 and S3. 

Treatment-dependent NAb development 

In agreement with previous studies [9,11,15,35], the frequency 
of NAbs was highest for IFNP-lb and lowest for i.m. IFNP-la, 
whereas high titers were most frequendy induced by s.c. IFNP-la, 
regardless of allele carriage. In the Swedish NAb registry the 
frequency of NAb development, independent of HLA allele group 
carriage, was 9.0% for i.m. IFNP-la, 33.6% for s.c. IFNP-la and 
49.3% for IFNP-lb. The frequency of biologically relevant titers 
was 3.9% for i.m. IFNP-la, 19.0% for s.c. IFNP-la and 15.9% for 
IFNP-lb (Table S4). 

HLA allele associations in the three different treatment 
groups 

Since NAb development is influenced by the IFNP preparation 
we analyzed association of HLA alleles within strata based on the 
IFNP preparation used. The AR of NAb development was also 
calculated for each allele group to assess how HLA allele group 
carriage influenced the risk of developing NAbs. 

After correction for multiple testing DRB1*15 carriage was 
associated with biologically relevant titers (OR 4.36; P c = 0.044; 
AR 5.8%) for patients receiving i.m. IFNP-la (Table 3). All 
analyzed HLA allele groups can be found in Table S5. 

For patients receiving s.c. IFNP- la an association was found for 
DRB1*15 with NAb positivity (OR 4.15; P,<0.001; AR 43.3%) 
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and with biologically relevant titers (OR 8.16; P„<0.001; AR 
27.7%) (Table 3). Associations were also found for carriage of 
DRB1*15 haplotype alleles C*07, B*07, DQB1*06 and DQA1*01 
(Table 3), which consequently disappeared upon DRB1*15 
stratification. In addition, DRB1*08 was found to be associated 
with NAb development (OR 3.42; P c = 0.036; AR 60.5%), 
however an association was not found within the smaller group 
of patients who developed biologically relevant titers (Table 3). 
Carriage of DQA1*05, which is commonly found on the same 
haplotype as DRB1*03, was associated with a decreased risk of 
developing biologically relevant titers (OR 0.29; P c = 0.032; AR 
8.4%) (Table 3). Carriage of DRB1*04 and DQB1*03 were only 
nominally associated with NAb development and could represent 
the same signal because of extensive LD between these allele 
groups. 

For patients receiving IFNP- lb, DRB1*04 carriage was 
associated with development of biologically relevant titers (OR 
3.53; P= 0.029; AR 28.2%) (Table 3). An association to HLA- 
DRB1 *04 with all NAb positive patients regardless of titer level 
could not be replicated in the Swedish cohort. Since an association 
has previously been seen between DRB1*04 subtypes and NAb 
development, we extended our analysis for biologically relevant 
titers to also include DRB1*04 subtypes (Table S6). The OR for 
NAb development in DRB1*04:01 carriers compared to non- 
carriers was 3.43 (95% CI 1.48-7.93; P= 0.0039), while no 
association was found for DRB1 *04:04 carriage (Table S6). 

Although the influence of allele group carriage was taken into 
account when estimating the risk of developing NAbs for each 
specific treatment the treatment adjusted AR for NAb develop- 
ment was constantly higher for the two s.c. IFNP preparations 
than for the i.m. IFNP preparation, regardless of whether the 
respective risk alleles DRB1*15 for IFNP- la and DRB1*04 for 
IFNP- lb were present or absent (Figure 1). For DRB1*15 carriers 
the AR of developing NAbs was higher if receiving s.c. IFNP- la 
(43% vs. 16%) and lower if receiving IFNP-lb (60% vs. 42%) 
compared to DRB1*15 non-carriers, while the opposite was 
observed for DRB1*04 carrier-ship (Figure 1, Table 3). Similar 
results were observed for biologically relevant titers. Although the 
risk of NAbs for i.m. IFNP- la was much lower it followed the 
pattern of s.c. IFNP- la (Figure 1, Table 3). 

Binding prediction of HLA alleles to IFN(3 peptides 

The preparation-specific association to different HLA alleles 
was further investigated by analyzing differences in predicted 
binding of HLA alleles to IFNP-la and IFNP-lb peptides. 
Especially the sites in IFNP that differ between the two IFNP 
molecules were investigated. The DRB1*04:01 allele, which was 
associated with an increased risk for NAbs in IFNP-lb treated 
patients, was found to bind weakly to IFNP-lb peptides containing 
Ser 17 but not to IFNP-la peptides containing Cys 17 . Similar 
distinctions between the two IFNP molecules could not be made 
for the alleles DRB1*04:04 and DRB1*15:01. Furthermore, 
DRB1*04:01 was predicted to bind with high affinity to peptides 
that included Asn 80 . However, this site is glycosylated in IFNP- la 
which could interfere with the binding to HLA alleles. Two IFNP- 
la peptides containing Meti were found to bind strongly to HLA- 
0*07:01. 

Discussion 

Preparation-related factors like dosage, injection site and 
frequency, aggregate content and genetic factors [9,16,20-22] 
may all contribute to IFNP immunogenicity. In this cohort we 
have investigated the role of HLA gene carriage for the risk of 



developing overall NAbs and titers of NAbs high enough to affect 
the biological response to IFNp. After stratification for the type of 
IFNP molecules (IFNP-la and IFNP-lb), we observed genetically 
distinct associations, DRB1*15 for IFNP-la and DRB1*04 for 
IFNP-lb. 

The risk for high NAb titers overall was 1.6-fold in patients 
carrying the MS risk haplotype DRB1*15, and for patients 
receiving IFNP- la the risk was increased up to 8-fold. The 
association between DRB1*15 and development of NAbs as well as 
high NAb titers is consistent with a previous study describing an 
increased proliferative response of CD4+ T cells against IFNP 
epitopes in DRB1*15 carriers compared to non-carriers [18]. 
CD4+ T cell activation is important for induction of immunolog- 
ical memory and persistent B-cell responses [36] and the 
interaction between T cells and B cells through HLA class II 
suggests that HLA might influence this type of response. 

There have been a few studies describing HLA associations with 
development of anti-drug antibodies against IFNP [20-22]. In the 
study by Hoffmann et al, which included an initial study cohort 
and a validation cohort, DRB1*04:01 and DRB1*04:08 were found 
to be associated with an increased risk of developing antibodies 
against IFNP in both groups [21]. These associations were further 
confirmed by Buck et al., using in an independent validation 
cohort [20]. Unlike these studies we found an association between 
NAb development and DRB1*15, whereas an association to 
DRB1*04 only was observed in patients treated with IFNP-lb. 
Factors that could contribute to the discrepancy in genetic 
association between the studies are the differences in allele 
frequencies between the Swedish and the German populations 
regarding DRB1*15 (15.6% in Sweden vs. 14.2% in Germany) and 
DRB1*04 (18.9% vs 13.2% in Germany, according to http:// 
www.allelefrequencies.net/default.asp; accessed December 13 in 
2013) as well as the different proportions of patients receiving 
IFNP-la (78% vs 60%) and IFNP-lb (22% vs 40%) in our cohort 
compared to the German cohort [21]. The fact that our cohort 
consisted of a higher proportion of IFNP- la users in which NAb 
development was associated to DRB1*15, and a lower proportion 
of IFNP-lb users in which NAb development was associated to 
DRB1*04, could in part explain why we only observed an 
association to DRB1*15 in the total cohort of this study. 
Additionally, in the previous studies, both binding and neutralizing 
anti-drug antibodies against IFNP were assessed together [20,21] 
whereas in the present study only patients who were positive for 
neutralizing anti-drug antibodies against IFNP were included. 
However, in contrast to the studies by Hoffmann et al. and Buck et 
al., we did not utilize a separate cohort for validation of our 
findings, which is a limitation of this study and thus these results 
need to be replicated in an independent cohort. The analysis of a 
possible association between DRB1 *04:08 and NAb development 
could not be performed due to the origin of the HLA data. The 
problematic assignment of alleles to patients with imputed 
genotypes could potentially influence our results. However, 92% 
(33 of 36 patients) of the known DRB1*04:08 positive individuals 
from another cohort [30] were imputed as DRB1*04:01 when 
investigating this problem, and carriage of both DRB1*()4:01 and 
DRB1*(>4:08 alleles has been shown to confer a higher risk of 
developing antibodies against IFNP [20,21]. 

This is the first report of distinct DRB1 allele associations for the 
two different types of recombinant IFNP molecules, DRB1*15 for 
IFNP-la and DRB1*04 for IFNP-lb, indicating an IFNP molecule 
dependent genetically determined risk to develop NAbs against 
IFNp. There are several possible explanations that could account 
for differences in immunogenic potential. Sequence alterations, 
different formulation contents with a potential to evoke the 
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Figure 1 . Absolute risk for development of NAbs and biologically relevant titers. Calculation of the absolute risk was used to estimate how 
DRB1*04 and DRB1*15 carriage impacts the risk for the outcomes of NAb positivity (bars on grey floor) and biologically relevant titers (bars on white 
floor) when adjusted for the baseline risk. Baseline risk was the frequency of development of NAb (dark grey bars, grey floor) and biologically relevant 
titers (dark grey bars, white floor) for each IFNp preparation in the Swedish NAb registry. Compared to the baseline risks, DRB1*15 carriage increased 
the absolute risk for both outcomes in s.c. IFNp-1a treated patients (grey bars) and to a lesser extent in patients receiving i.m. IFNp-1a (light grey 
bars), whereas DRB1*15 carriage lowered the risk for both outcomes in patients receiving IFN(3-1b (white bars). In DRB1*04 carriers the reversed 
relationship was observed, with an increased absolute risk for both outcomes in IFN(3-1b treated patients. Moreover, the absolute risk for NAb 
positivity and biologically relevant titers was constantly lower for i.m. IFNf3-1a than for the two s.c. IFN(3 preparations, regardless of whether patients 
were positive or negative for DRB1*04 and DRB1*15. 
doi:1 0.1 371 /journal.pone.0090479.g001 



immunological response, and different propensities to form 
aggregates due to the glycosylation status are all factors that 
might have a higher impact on NAb development. The different 
HLA associations found here might be explained by different 
binding affinities to the two types of IFNP molecules to different 
HLA alleles. For instance, the increased risk for NAbs against 
IFNP- lb in DRB1*04:01 carriers compared to non-carriers might 
to some extent be explained by the differential binding potential of 
this allele to IFNP- lb compared to IFNP- la. 

Binding prediction is based on the primary IFNP sequence, but 
it is also possible that post-translational modifications such as 
glycosylation might interfere with peptide binding and presenta- 
tion in vivo. Naturally processed peptides containing Asn-linked 
glycosylations can be presented on HLA class II molecules and 
these glycans can also alter the processing and presentation of 
peptides [37,38]. We only observed an association with NAb to 
DRB1*04:01 in patients treated with IFNP- lb, the non-glycosy- 
lated preparation. One could thus hypothesize that DRB1*04:01 
carriage might lead to more efficient presentation of non- 
glycosylated peptides. This might be one explanation for the 
different HLA association to different IFNP molecules. 

Although HLA alleles influence the risk of NAbs within the 
treatment groups it is clear that NAb development in MS patients 
treated with IFNP- la is also influenced by several additional 
factors. Thus, IFNP- la injected intramuscularly is still the least 
immunogenic way of receiving IFNP even for DRB1*15 carriers 
and none of the analyzed HLA allele groups lowered the risk of 
subcutaneously injected IFNP to the level of intramuscularly 
injected IFNP- la. The skin harbors a high frequency of dendritic 



cells efficient in presenting antigens [39], easily provoked by 
repeated injections of IFNP and potentially explaining higher 
immunogenicity of the subcutaneously administrated IFNP 
products compared to the intramuscularly administered product. 

Approximately 20% of IFNP-treated patients develop NAb titers 
high enough to block treatment bio-efficacy. Since up to 60% of MS 
patients carry DRB1*15 and some 30% carry DRB1*04, and since 
the risk of high titer NAbs is increased several times (ORs 3-4) in 
carriers of respective risk alleles, it may well be relevant to base 
treatment decisions on HLA genotype. In fact, our empirical data 
suggest that DRB1*15 carriers would have a 35% risk of losing 
efficacy of s.c. IFNP- la but less than 10% on IFNP- lb (provided 
they did not carry DRB1*04 on their other chromosome). The 
reverse risks would be seen for a DRBl*15-ntga.txvt carrier, who 
carries DRB1*04, on IFNP- lb. Although this study is based on over 
900 patients the observed ORs show, however, large confidence 
intervals, making such clinical procedures somewhat premature. 
Extended studies on even larger patient cohorts are therefore 
warranted. Our findings may, though, have some bearing on other 
biological treatments for other diseases where blocking or neutral- 
izing antibodies pose a problem. 

Conclusion 

HLA class II genotype influences the risk of MS patients to 
develop NAbs against IFNP-la and IFNP-lb. Further studies of 
the genetic risks for induction of anti-drug antibodies can provide 
important knowledge on possible immunological mechanisms 
triggered by regularly injected biopharmaceuticals. 
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